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We extend our investigation of heavy-light meson-meson interactions to a system consisting of a heavy-light me- 
son and the corresponding antiparticle. An effective potential is obtained from meson-antimeson Green-functions 
computed in a quenched simulation with staggered fermions. Comparisons with a simulation using an 0(a 2 ) 
tree-level and tadpole improved gauge action and a full QCD simulation show that lattice discretization errors 
and dynamical quarks have no drastic influence. Calculations from inverse scattering theory propose a similar 
shape for KK potentials. 



Motivated by the success of QCD in high en- 
ergy scattering, today one aim is to calculate 
hadron-hadron interactions directly from the field 
equations of QCD. In the low energy regime of 
QCD nonperturbative tools have to be used. This 
leads us to study the forces in systems of two 
hadrons on the lattice. 

Previous lattice calculations have revealed a 
short range attractive potential between two iden- 
tical heavy- light color singlets In this paper 
we present a lattice-QCD study of a heavy-light 
meson-antimeson system MM and compare the 
resulting residual interaction with KK potentials 
obtained from inverse scattering theory. 

We define the one-meson field as a product of 
staggered Grassmann fields x and X with a heavy 
and a light external flavor h and I, respectively, 



in the time correlation matrix 

CW,(Mo) = (4(t)^(to))-<4(*))(^'(to)),(3) 

where ( ) denotes the gauge field configuration 
average. On the hadronic level C is a two-point 
correlator of a composite local operator describ- 
ing a molecule-like structure. Working out the 
contractions between the Grassmann fields we ob- 
tain 



Cp\t, to) 



xt , 



^i G x+rt,x+rt a G x+rt,x+rt )) 

+ .... 



(4) 



4>s{t) = Xh{xt)xi{xt) . 



(1) 



The meson-antimeson fields with relative distance 
r = y — x are then constructed by 



Since the heavy valence quarks are fixed in space 
the relative distance between the mesons is the 
same at the initial and final time of the propaga- 
tion, r ' = f. The heavy-quark propagator is 



(2) 



C {h) = 
xt,xt n 



1 



[T S xi] k Y\_ ^=(ij'a),(i=4 ! (5) 



x y 



The dynamics of the MM system is contained 
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where the phase factors T S4 = (_l)(*i-Hca+x a )/o 
in the Kogut-Susskind formulation correspond to 
the Dirac matrices and k = (t — to) /a. In our 
calculations we set m^a = 0.5. The propagator 
of the light quark is obtained from inverting the 
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staggered fermion matrix with a random source 
estimator. A standard conjugate gradient algo- 
rithm is used. 

The effective ground-state energy of the MM 
system W(r) can be extracted from the large eu- 
clidean time behavior of C following quantum- 
mechanical reasoning for composite particles, 

The residual meson- antimeson potential is 
V{r) =W(r) -2m, (7) 

where 2m is the mass of two free (anti)mesons. 

The simulations were performed on a periodic 
8 3 x 16 lattice. Each potential is the result of 
a measurement on 100 independent gauge field 
configurations separated by 200 sweeps. The in- 
version of the fermion matrix was performed with 
32 random sources. 

Each contribution to the correlator in (Q) com- 
prises the exchange of gluons. As a first step only 
the direct term C^ A ' was computed, the com- 
putation of the full correlator being the task of 
a subsequent work. The effective energy W{r) 
was then extracted from the correlator by 
a four parameter (Levenberg-Marquardt) fit with 
the function 

C [ r A) = B(r) cosh [W(r) (t - 8a)} 

+ (-iy/ a B(r)cosh[W(r)(t-8a)]. (8) 

The second term alternating in sign is a peculiar- 
ity of the staggered scheme. To be numerically 
consistent with W(r) at large distances, the mass 
2m of two noninteracting mesons was extracted 
from the square of the one-meson two-point func- 
tion. 

In order to investigate the influence of dis- 
cretization errors and sea quarks a simulation us- 
ing an improved gauge action and a full QCD sim- 
ulation were also performed. The former includes 
planar six-link plaquettes in addition to the ele- 
mentary plaquettes and is tadpole improved B]. 
The resulting potentials are displayed in Fig.^l]. 
All simulations show consistent results. The po- 
tential turns out to be short ranged and attrac- 
tive. At distances r/a > 1 essentially no interac- 
tion energy can be resolved. 
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Figure 1. Heavy- light meson-antimeson poten- 
tials V{r)a for quark mass m/cs = 0.1, from a 
quenched simulation using the Wilson gauge ac- 
tion with = 5.6, from an improved simulation 
with P p i — 7.0 corresponding to the same lattice 
constant a w 0.2 fm, and from a full QCD simula- 
tion with number of flavors n/ = 3 and /3 = 5.2. 



Data from an independent simulation for a dif- 
ferent light-quark mass mfa — 0.05 were used to 
make a linear extrapolation to light quark mass 
zero. To estimate the strength of the interac- 
tion the data were matched to a Gaussian po- 
tential. The latter was then used to obtain scat- 
tering phase shifts from a Schrodinger equation 
with three arbitrary values of the meson mass. 
The results are shown in the left plots of Fig. |[ 

An effective range expansion with two sets of 
parameters was applied to experimental data of 
KK s-wave scattering || . Using quantum inver- 
sion these analytic phase functions can be trans- 
formed into local potentials . The input phase 
shifts and the resulting potentials are shown in 
the right plots of Fig. [| We notice a good qualita- 
tive agreement between lattice and inversion po- 
tentials for short distances. The repulsive hump 
could not be resolved by our simulation. This 
feature, as it turns out, is decisive for the phase 
shifts and explains the lack of qualitative agree- 
ment between Sq. 

The next step in our program is to compute 
the full correlator C r . A new simulation with an 
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Figure 2. Potentials and phase shifts for the MM system from the lattice (left plots) and for the KK 
system from inverse scattering (right plots). 



0(a 2 ) tree- level and tadpole improved gauge ac- 
tion and Naik-type improved staggered fcrmionic 
action is in progress. An independent lattice 
investigation with 0(a 2 )-improved lattice action 
using Wilson fermions is presented in these pro- 
ceedings. 
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